Some aspects of the critical behavior of the 
Two-Neighbor Stochastic Cellular Automata 



i/^ ' A.Yu. Tretyakov 

School of Information Sciences, Tohoku University, 
Aramaki, Aoba-ku, Sendai, 980-77, JAPAN 

Oh! 

$ : N. Inui 

! Himeji Institute of Technology, 

2167, Shosha, Himeji, 671-22, JAPAN 

> ■ M. Katori and H. Tsukahara 

Department of Physics, Faculty of Science and Engineering, Chuo University, 

§n : Kasuga, Bunkyo-ku, Tokyo 112, JAPAN 

O ■ 
in : 

^ ■ February 6, 2008 

-i— » ■ 
ctf ■ 

a : 

i 1 Abstract 

■ 

Using Pade approximations and Monte Carlo simulations, we study the phase diagram 
of the Two-Neighbor Stochastic Cellular Automata, which have two parameters pi and 

•i-H . 

^ ■ p2 and include the mixed site-bond directed percolation (DP) as a special case. The 

phase transition line pi = picfe) has endpoints at (pi,P2) = (1/2, 1) and at (0.8092, 0). 
The former point (1/2,1) is a special point at which Compact DP transition occurs and 
its critical exponents are known exactly. Results of time-dependent simulation show that 
in the whole range of parameters, excluding this point (1/2,1), the system belongs to the 
DP universality class. It is first shown that the shape of the phase transition line near 
this special point has, asymptotically, a parabolic shape, i.e., Pi c (p2) — 1/2 ~ (1 — P2) 6 
with 9 = 1/2 for < 1 — pi <S 1. We use the Monte Carlo data to assess the accuracy of 
rigorous bounds for the line recently reported by Liggett and by Katori and Tsukahara. 
It is also shown that outside the vicinity of the special point (1/2,1), the curve is well 
approximated by an interpolation formula, similar to the one proposed by Yanuka and 
Englman. 
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1 Introduction 



The two-neighbor stochastic cellular automata (SCA), introduced by Domany and Kinzel, can 
be regarded as a very general model of a spread of influence in 1+1 dimensions. A site on layer 
% + 1 is occupied with probability pi if one and only one of its nearest neighbors on the previous 
layer is occupied, it is occupied with probability p 2 if both of the neighbors are occupied, and 
it remains vacant if both neighbors are vacant (Fig. p]). Often the system is regarded as a 
process taking place in one-dimension, with the layer number % replaced by discrete time t. 

As one can easily see, the case of p% = p 2 = p is equivalent to the site directed percolation (DP) 
with site occupation probability p. It is easy to show that (p± = p, p 2 = 2p — p 2 ) corresponds 
to bond DP with the fraction of open bonds p. The mixed site-bond DP with the fraction 
of open sites a and the fraction of open bonds ft corresponds to the two-neighbor SCA with 
pi = aft and p 2 = a{2ft - ft 2 ) |, |, |. 

In the present work, by using Monte Carlo simulations, we evaluate the phase transition line 
which is characterized as follows. We consider the process starting from a single occupied 
particle and let P(t) be the probability that at least one particle survives at time t. In the 
subcritical phase p\ < Pi c (p2), this survival probability decreases exponentially to zero, while 
in the supercritical phase pi > pi c (p2), it converges to a positive value: P^ = lim^oo P(t) > 0. 
On the phase transition line p\ = Pi c (p2) ; P{t) shows a power-law decay and critical phenomena 
are observed. The ultimate survival probability in the supercritical phase behaves as follows 
in the vicinity of the phase transition line: P 00 (p 1 ,p 2 ) ~ (Pi — Picfe))^ as pi — > Pi c (p2) for 
each p 2 . It is conjectured for p 2 < 1 that the critical exponent ft is universal, i. e., independent 
of p 2 and is the same as the pure DP, ft ~ 0.276. 

The two-neighbor SCA with p 2 = 1 is equivalent to the reaction-limited case of the so-called 
dimer-dimer model Q, this case (sometimes referred to as Compact Directed Percolation) is 
understood rather well [3], In particular, it is easy to see P oc (p 1 , 1) = 1 — (1 — pi) 2 /Pi — 
8(pi — 1/2) for pi ~ 1/2, which means p lc = 1/2 and ft(p 2 = 1) = 1 [[J. 



The point (pi,p 2 ) = (1/2, 1) is an endpoint of the critical line pi = pi c (p2) i n the phase diagram 
of the two-neighbor SCA. The universality for p 2 < 1 and the above mentioned exact result for 
P2 — 1 imply that we will observe the following phenomenon called a crossover when p 2 — 1 @ . 
There is a function e{p 2 ) > such that if < p\ — Pxdp'i) "C s{p2) then P^ ~ (pi — picfe))' 3 
with (3 ~ 0.276, while if e{p 2 ) < P\ — Picfe) <^ 1 then ~ (j>i — 1/2) and that e(p2) ^ as 
P2 -> 1. 

In order to study the universality and the crossover, it is necessary to make precise evaluation 
of the phase transition line. In particular, we need to know the asymptotic behavior of the 
phase transition line near the special point (pi,p 2 ) = (1/2,1). By using series expansion 
technique and comparing the results with the data obtained by Monte Carlo we conclude that 
the asymptotic shape of the phase transition line near (1/2, 1) is parabolic. In other words, 
Pi c (p 2 ) - 1/2 ~ (1 - p 2 ) e with 6 = 1/2 for < 1 - p 2 < 1. We use the Monte Carlo data to 
assess the accuracy of rigorous bounds recently obtained theoretically || [?]]. 

2 Monte Carlo simulations 



Two points of the phase diagram, corresponding to pure site and bond DP, can be estimated 
from the corresponding threshold values as (0.705485,0.705485) and (0.644701,0.873762), using 
the data taken from M. The sector of the phase diagram between them corresponds to the 



mixed site-bond DP and can be obtained using the data given in [12]. Therefore, to obtain a 
complete phase diagram only areas with 0.5 < p\ < 0.644701 and p\ > 0.705485 need to be 
covered. 



Similarly to [12[], to determine the location of the phase-transition line, we used the time 
dependent simulation technique and analyzed the time dependence of the local slopes for the 
mean number of particles n, for the survival probability P and for the mean square radius 
R. Simulations were conducted starting from a single occupied particle, taking 1000 steps and 
averaging over 10 5 configurations. Some of the local slopes curves (for the mean number of 



particles) are given in Fig. [|. Our results agree with the assumption that the two-neighbor 
SCA belong to the DP universality class for p 2 < 1 |J, with the dynamical critical exponents 
defined as 

P(t) oc t- s 

n{t) oc t" (1) 
R 2 {t) oc t z 

falling within intervals 5 = 0.162 ± 0.020, rj = 0.308 ± 0.015 and z = 1.263 ± 0.040 from 
the DP values (we take the DP values of the dynamical critical exponents from |§). On the 
other hand, the local slopes approach the asymptotic values, corresponding to the values of 
the critical exponents, in a different way depending on the value of p2- While for the part 
of the phase diagram far enough from the special point (1/2,1) the local slopes stabilize at 
the asymptotic value relatively fast (Fig. |^a), for the points close to (1/2,1), in the cases we 
studied, the number of steps taken was not enough to reach the asymptotic value, although the 
approach to DP values was still quite obvious (Fig. §b). For (pi,p2) = (1/2, 1) the exponents 
agreed with the theoretical predictions of 5 = 1/2, 77 = and 2 = 1 |. 

The Monte Carlo results for the phase transition line are presented in Tab. 3. In Fig. § 
we compare them with the rigorous lower bound given by Katori and Tsukahara as P2 = 
(1 — 2pf)/(l — pf) and with the rigorous lower (p 2 = 2(1 — pi)) and upper (p 2 = ^Pi ~ 4p?) 
bounds given by Liggett 0. Also, we present the rigorous lower bound obtained by averaging 
the number of particles after 12 steps (obtained by exact enumeration, starting from a single 
particle) and taking pi and P2 values leading to the average equal to 1. The upper bound given 
by Liggett comes very near to the phase transition line in the vicinity of (1/2,1). It should 
be remarked that all these bounds can only be proved to be valid in the so-called attractive 
region, with pi < P2 [0]. However, we see that the Katori and Tsukahara's lower bound is not 
only valid over the whole range of P2 but, indeed, follows the shape of the Monte Carlo data 
rather closely. 

Fig. 3 should be compared with the phase diagram shown by Kinzel [Q], which was studied by 



the transfer-matrix scaling method. Kinzel showed two possibilities for the behavior of the 
phase transition line near p 2 = 0: the line has an endpoint on the pi-axis at (pio,0) with 
Pio < 1 with pio < 1, or the line goes to the corner point (1,0) in the phase diagram. The 
former case is concluded from our results with p w = 0.8092 ± 0.0004. 



An interpolation formula, proposed by Yanuka and Englman |TT[ for the ordinary percolation, 
was shown to perform rather well in the case of the mixed DP [0. Expanding it to the general 
case of the two- neighbor SCA leads to 

p 2 = Pi(2-( Pl M) 1/(1 - A) )- (2) 

where A = logpf/logp|, with p\ and p\ as ordinates of points on the phase diagram, corre- 
sponding to the pure site and pure bond DP, respectively. Although we can not present any 
theoretical explanation at this moment, one can see that (0) perform rather well even outside 
the mixed site-bond percolation sector, for which it has been derived. Marked discrepancy 
begins at p 2 > 0.9, as the phase diagram approaches the special point (1/2,1). 

3 Series expansion 

While an attempt to directly apply Pade analysis to the series expansion for the survival 
probability in p\ and P2 around (0,0) fails because all coefficients of the expansion change 
when the order is increased, changing coordinates to a and q according to 

Pi = 1 - aq (3) 
p 2 = 1-q (4) 

leads to a series with converging coefficients. Computer time required to obtain, using exact 
enumeration, an expansion to a certain order grows very fast with the order of the expansion 
because the model includes 2 parameters. In the current work we try to take advantage of the 
possibility to work analytically with a relatively low order expansion. 



Regarding the series expansion for the logarithm of the survival probability in a and q as a 
polynomial in q we take a [3,3] Pade approximant for its derivative in q and expect that the 
phase transition line is given by a root of the denominator (for a description of the DLog Pade 
method we are using, for example, fL3|). Solving the corresponding 3-rd order equation exactly 
and discarding unphysical roots we get a rather complicated expression for the phase diagram 
in the form of q given as a function of a. Assuming that we are in the vicinity of the special 
point, where, as one can easily see, a goes to infinity, we take an expansion in 1/a 1 / 2 , the two 
leading terms are obtained as 

where the discarded terms are of the order 0(a 5//2 ). As one can see in Fig |], ([5]) approximates 
the exact solution surprisingly well. Taking @ and returning to coordinates p± and p 2 we have 
for the phase diagram in the vicinity of (1/2,1) 

p 2 = -3 + 20 Pl - 32pj + 16pl (6) 

Again, taking a Taylor expansion in pi — 1/2 we have 

P2 = 1 - 8(px - 1/2) 2 , (7) 

discarding terms of the order 0((pi — 1/2) 3 ). Surprisingly enough, comparison with the Monte 
Carlo data (Fig. ^) shows that both and (0) work reasonably well in the vicinity of (1/2,1), 
with (^), as it should be expected, working over a wider range of p\ values. Moreover, the 
parabolic asymptotic behavior, predicted by (|7|), is in agreement with the Ligget's upper bound 
and with the Monte Carlo data, as it is shown in Fig. [5| 

4 Conclusions and Future Problems 



By using Pade approximations and Monte Carlo simulations, we evaluate the phase transition 
line P i = Pic{p2) for the two-neighbor SCA. The line has endpoints at (pi,P2) = (1/2, 1) and at 
(0.8092,0). Using the time dependent simulation technique Hl2 |, we confirm that the present 



SCA belong to the DP universality class in the whole range of parameters, excluding the point 
(1/2,1). The special point (1/2,1) is the transition point of the reaction-limited case of the 
dimer-dimer model and different critical phenomena are observed [[J. The present work 
show that the shape of the phase transition line near this special point has, asymptotically, a 
parabolic shape, and we conclude that 

PiM - ^ ~ C(l - p 2 ) e with 6 = ^ and C ~ 0.436 (8) 

for < 1 — P2 <C 1 . As discussed in Section 1 , this numerical result is a first step in the further 
study on the crossover which is expected to occur in the vicinity of (^1,^2) = (1/2, 1) ||. 

Recently rigorous lower and upper bounds of the phase transition line were given for the 
attractive region p\ <p2 H 0- In this paper we present other rigorous lower bound obtained 
by exactly averaging the number of particles after 12 steps, whose validity can also be proved 
if and only if p 1 <P2- We access the accuracy of these rigorous bounds using the Monte Carlo 
date. The upper bound given by Liggett is excellent in the vicinity of (1/2,1). It is shown 
that the lower bound by Katori and Tsukahara seems to be valid also in the non-attractive 
region p\ > p2 and follows the shape of the Monte Carlo date rather closely down to P2 = 0. 
Justification of the "lower bound" for pi > p 2 is a challenging problem, since rigorous things 



proved for the non- attractive region are still limited flT4|. 



We also show in this paper that far enough from the special point (1/2,1) the phase transition 



line is well approximated by an interpolation formula [pL2|| , similar to the one proposed by 
Yanuka and Englman for ordinary percolation [|TT . 
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Tables 



1. Monte Carlo results for the phase diagram. 

Values marked by * were obtained with uncertainty ±0.0002, ' corresponds to uncertainty 
±0.0004, and * to ±0.0008. 

a) the vicinity of p 2 = 1 . 

b) the vicinity of p 2 = 0. 

Figure captions 

1. The Two-Neighbor Stochastic Cellular Automata 

2. Local slopes for the mean number of particles r = \og(n(t)/n(t/5))/ log(5) plotted versus 
the inverse of the time t. 

Horizontal solid line corresponds to the DP value of the critical exponent r\ = 0.308. 

a) p 2 = 0.62, 1-px = 0.9175, 2-pi = 0.9171, 3-pi = 0.9167 (critical point), 4-pi = 0.9163, 
5- pi = 0.9159. 

b) p 2 = 0.52, 1- Pl = 0.9983, 2- Pl = 0.9981, 3- p x = 0.9979 (critical point), 4- Pl = 
0.9977, 5- pi = 0.9975. 

c) pi = 1/2, p 2 = 1 

3. Phase diagram for the Two-Neighbor SCA. Monte Carlo results are represented by points. 
1,2- lower and upper bounds by Liggett 0. 

3- lower bound by Katori and Tsukahara 

4- Pade result. 

5- Tailor expansion for the Pade result. 



6- Interpolation formula [11, 12 



7- lower bound obtained by averaging the number of particles after 12 steps (obtained by 
exact enumeration, starting from a single particle) and taking pi and p 2 values leading 
to the average equal to 1. 



4. Phase diagram by Pade approximants in coordinates q and a, defined by p\ = 1 — q and 
p 2 = 1 — aq. The lower line corresponds to the Taylor expansion, given by q = (| — \-)/a. 

5. Phase diagram for the Two-Neighbor SCA in double logarithmic coordinates. Solid line 
shows the best fit for q = 1 — p 1 = k{pi — 0.5) 2 over the first 3 points, which corresponds 
to k = 5.25 . . .. Taylor expansions for the Pade result (k = 8) and for the Liggett's upper 
boundary {k = 4) are also shown. 



i n 



Pi P2 

0.50 1.0000 

0.52 0.9979* 

0.54 0.9917* 

0.56 0.9809* 

0.58 0.9653* 

0.60 0.9443* 

0.62 0.9167 t 



Pi P2 

0.71 0.6877 t 

0.73 0.6010+ 

0.75 0.4933+ 

0.77 0.3590+ 

0.79 0.1916+ 

0.8092+ 0.0000 



1 1 



